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Abstract 


The  report  assesses  the  capability  of  Computational  Fluid  Dynamics  to  predict  the  aerother- 
modynamic  loads  on  simplified  geometries  in  hypersonic  flight.  The  specific  configurations 
selected  are  the  double  cone  and  hollow  cylinder  flare.  Experiments  were  previously  con¬ 
ducted  at  the  Calspan  University  of  Buffalo  Research  Center  (CUBRC)  at  stagnation  en¬ 
thalpies  from  5  MJ/kg  to  22  MJ/kg  in  air  and  Mach  numbers  from  10.9  to  13.2.  Experimental 
diagnostics  include  surface  pressure  and  heat  transfer.  Laminar  compressible  Navier-Stokes 
simulations  were  performed  using  both  perfect  gas  and  non-equilibrium  Park  I  thermochem¬ 
istry  models.  Detailed  grid  refinement  studies  were  performed  to  establish  the  accuracy  of 
the  simulations.  The  principal  conclusions  of  the  research  are 

•  The  Park  I  model  accurately  predicts  the  peak  heat  flux  max  Qu  in  nearly  all  cases 

•  The  Park  I  model  accurately  predicts  the  separation  length  Lsep  and  peak  pressure 
max  p  in  approximately  one-half  of  the  cases 

•  The  perfect  gas  model  accurately  predicts  the  peak  heat  flux  max  Qw  in  all  four  hollow 
cylinder  flare  cases,  but  in  none  of  the  double  cone  cases 

•  The  perfect  gas  model  accurately  predicts  the  separation  length  Lsop  in  one  case  and 
peak  pressure  max  p  in  none  of  the  cases 


The  significance  of  the  results  is: 

•  The  Park  I  model  is  a  reliable  model  for  prediction  of  peak  heat  flux 

•  The  Park  I  model  is  not  a  reliable  model  for  prediction  of  peak  surface  pressure  or 
separation  length  (ie.,  surface  pressure  distribution) 

•  The  perfect  gas  model  is  not  a  reliable  model  for  prediction  of  peak  surface  pressure  or 
separation  length  (i.e.,  surface  pressure  distribution)  except  at  the  lowest  stagnation 
enthalpy 

•  The  perfect  gas  model  is  not  a  reliable  model  for  prediction  of  peak  heat  flux  except  at 
the  lowest  stagnation  enthalpy;  it  is  accurate  for  hollow  cylinder  flare  at  all  stagnation 
enthalpies,  but  not  for  the  double  cone 

The  significance  is  limited  to  the  types  of  shock  wave  laminar  boundary  layer  interactions 
considered  in  this  study. 
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1  Introduction 


The  design  of  effective  high  speed  ( e.g .,  hypersonic)  weapons  systems  requires  accurate 
modeling  of  the  influence  of  the  flow  on  the  vehicle  performance.  In  particular,  the  interaction 
of  shock  waves  and  boundary  layers  leads  to  a  variety  of  critical  flow  phenomena  including 
high  heat  transfer  at  reattaclunent,  flow  separation  resulting  in  significant  changes  to  the 
surface  pressure  (and  lienee  forces  and  moments  on  the  vehicle),  and  unsteadiness  causing 
dynamic  aerothermodynamic  loads.  Thus,  accurate  modeling  of  shock  wave  boundary  layer 
interaction  is  essential. 

Shock  wave  boundary  layer  interaction  has  been  studied  since  the  1940s.  One  of  the  earliest 
experiments  was  performed  by  Liepmann[l]  examining  the  interactions  of  a  shock  wave  on 
a  twelve  percent  circular  arc  airfoil  at  Mach  0.8  for  both  laminar  and  turbulent  boundary 
layers.  Schlieren  imaging  (Fig.  1)  and  surface  pressure  measurements  indicated  a  significant 
effect  of  the  boundary  layer  state  (i.e.,  laminar  vs  turbulent). 


Figure  1:  Shock  wave  boundary  layer  interaction  [1]  at  Mach  0.863 


Beginning  in  the  late  1960s,  Computational  Fluid  Dynamics  (CFD)  emerged  as  a  separate 
branch  of  fluid  mechanics  research.  Early  on  significant  effort  was  focused  on  simulation 
of  high  speed  flows,  in  particular,  shock  wave  boundary  layer  interaction  (e.g.,  MacCor- 
mack  [2],  Shang  et  al  [3],  Hung  and  MacCormack  [4]).  Beginning  in  the  early  1990s,  CFD 
modeling  of  shock  wave  boundary  layer  interactions  had  progressed  to  a  degree  sufficient 
to  warrant  detailed  assessment  by  NATO  AGARI).  An  early  assessment  was  performed  by 
AGARD  Working  Group  18  (1992-1997)  focusing  on  shock  wave  boundary  layer  interaction 
at  supersonic  speeds.  Laminar  and  turbulent  boundary  layers  were  considered.  Only  RANS 
models  were  considered  in  the  latter  case.  The  study  highlighted  the  failure  of  RANS  model 
predictions  of  skin  friction  and  heat  transfer  (Knight  and  Degrez  [5]). 

AGARD  Working  Group  10  (1998-2003)  continued  evaluation  of  supersonic  shock  wave 
turbulent  boundary  layer  interactions.  Modeling  was  extended  to  include  LES  and  DNS 
(for  nominal  two-dimensional,  low  Reynolds  number  interactions)  and  RANS  (for  three- 
dimensional  interactions).  The  former  showed  significant  progress  in  predictive  capability; 
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however,  no  comparison  with  surface  heat  transfer  was  performed  and  the  Reynolds  numbers 
considered  were  not  representative  of  conditions  for  modern  air  vehicles.  The  RANS  models 
failed  to  accurately  predict  heat  transfer  (Knight  et  al  [6]). 

RTO  AVT  Task  Group  136  (2004-2008)  extended  the  focus  to  hypersonic  laminar  shock 
wave  boundary  layer  interactions  in  non-equilibrium  flows.  Two  configurations  were  con¬ 
sidered  -  a  double  cone  model  at  Mach  11.6  in  N2  at  stagnation  enthalpies  of  5.3  MJ/kg 
and  9.7  MJ/kg,  and  a  cylinder  in  air  at  Mach  9  at  stagnation  enthalpies  of  13.5  MJ  and 
22.4  MJ/kg.  CFD  capability  for  prediction  of  surface  heat  transfer  and  pressure  for  the 
double  cone  was  mixed;  poor  agreement  with  surface  heat  transfer  011  the  cylinder  at  the 
highest  enthalpy  was  observed  (Knight  et  al  [7]). 

STO  AVT  Task  Group  205  (2011-2014)  focused  on  hypersonic  laminar  shock  boundary  layer 
interaction  and  hypersonic  boundary  layer  transition  caused  by  a  single  protuberance.  CFD 
simulations  of  a  nominally  double  wedge  configuration  identified  the  limitations  of  the  ex¬ 
periment.  In  particular,  the  experimental  surface  heat  transfer  did  not  represent  the  fully 
steady  state  flow,  and  the  experimental  flowfield  was  evidently  three-dimensional.  (Knight 
and  Chazot  [8]). 

In  summary,  the  results  of  several  previous  assessments  of  CFD  capability  for  modeling 
hypersonic  shock  wave  boundary  layer  interaction  (in  particular,  for  laminar  boundary  layers 
at  high  stagnation  enthalpies)  is  inconclusive. 

The  principal  objective  of  this  project  was  to  assess  CFD  modeling  capability  for  hypersonic 
shock  wave  laminar  boundary  layer  interactions.  Two  experimental  configurations  were 
selected  -  the  double  cone  and  hollow  cylinder  flare  models  tested  at  the  Calspan  University  of 
Buffalo  Research  Center  (CUBRC)  at  stagnation  enthalpies  from  5  M  J /kg  to  22  MJ  /kg  in  air 
and  Mach  numbers  from  10.9  to  13.2.  Experimental  diagnostics  include  surface  pressure  and 
heat  transfer.  Laminar  compressible  Navier-Stokes  simulations  were  performed  using  both 
perfect  gas  and  non-equilibrium  Park  I  thermochemistry  models.  Detailed  grid  refinement 
studies  were  performed  to  establish  the  accuracy  of  the  simulations. 


2  Description  of  Experiments 

An  extensive  set  of  experiments  were  performed  at  Calspan  University  of  Buffalo  Research 
Center  (CUBRC)  under  the  direction  of  Dr.  Michael  Holden  to  examine  the  effects  of  thermo¬ 
chemistry  in  shock  wave  laminar  boundary  layer  interactions.  Two  separate  configurations 
were  considered  -  the  double  cone  (Fig.  2)  and  the  hollow  cylinder  flare  (Fig.  3).  The  exper¬ 
iments  were  performed  in  the  LENS  XX  facility.  A  description  of  the  facility  is  presented  in 
Dufrene  et  al  [9,  10].  The  inflow  gas  was  air  in  full  chemical  and  thermochemical  equilibrium 
with  mass  fractions  of  0.765  and  0.235  for  N2  and  O2,  respectively.  The  model  surface  was 
isothermal  at  300  K.  The  inflow  conditions  for  the  double  cone  and  hollow  cylinder  flare 
experiments  are  listed  in  Tables  1  and  2,  respectively. 
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(a)  Photograph 


Figure  2:  Double  Cone  (dimensions  in  inches  [mm]) 


iw  pn.«] 


(a)  Photograph 


(h)  Drawing 


Figure  3:  Hollow  Cylinder  Flare  (dimensions  in  inches  [mm]) 


Experimental  diagnostics  included  surface  lieat  transfer  and  pressure.  A  total  of  24  pressure 
transducers  and  50  heat  flux  gauges  were  installed  in  the  double  cone  model,  and  19  pressure 
transducers  and  52  heat  flux  gauges  in  the  hollow  cylinder  flare  model.  The  pressure  trans¬ 
ducer  is  (3  mm  diameter.  The  heat  flux  gauge  is  a  strip  whose  dimension  in  the  streamwise 
direction  is  less  than  C.l  mm.  Approximately  7%  of  the  pressure  transducer  data  and  4%  to 
12%  of  "che  heat  transfer  data  was  not  used.  The  experimental  uncertainty  for  heat  flux  and 
pressure  measurements  is  ±10%  [11]. 
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Table  1:  Flow  Conditions  for  Double  Cone  Experiments 


Run 

No. 

Total 

Enthalpy 

(MJ/kg) 

Mach 

Number 

Pitot 

Pressure 

(kPa) 

Unit 

Reynolds 

Number 

/106(m-1) 

Velocity 

(km/s) 

Density 

(g/m3) 

Temperature 

(K) 

1 

5.44 

12.2 

5.1 

0.14 

3.246 

0.499 

175 

2 

9.65 

10.90 

17.5 

0.19 

4.303 

0.984 

389 

6 

15.23 

11.46 

59.0 

0.39 

5.466 

2.045 

573 

4 

21.77 

12.82 

39.5 

0.20 

6.497 

0.964 

652 

NOTE:  “Run  No. 

”  corresponds  to  the  experimental 

run  identification  number 

Table  2 

:  Flow  Conditions  for  Hollow  Cylinder  Flare  Experiments 

Run 

No. 

Total 

Enthalpy 

(MJ/kg) 

Mach 

Number 

Pitot 

Pressure 

(kPa) 

Unit 

Reynolds 
(/106  m-1) 

Velocity 

(km/s) 

Density 

(gm/m3) 

Temperature 

(K) 

1 

5.07 

11.3 

5.9 

0.15 

3.123 

0.634 

189 

2 

10.43 

12.6 

9.7 

0.12 

4.497 

0.499 

318 

4 

15.54 

11.5 

64.0 

0.42 

5.470 

2.216 

569 

5 

21.85 

13.2 

39.0 

0.20 

6.515 

0.947 

618 

NOTE: 

“Run  No.” 

corresponds  to  the  experimental  ] 

run  identification  number 

In  comparing  the  computed  heat  transfer  and  pressure  with  experiment,  it  is  important  to 
relate  the  computational  grid  spacing  along  the  surface  to  the  sensor  size.  The  streamwise 
length  of  the  heat  sensor  is  less  than  100/zm,  and  the  diameter  of  the  pressure  sensor  is  6  min. 
The  grid  resolution  of  the  simulations  in  the  vicinity  of  the  peak  pressure  is  typically  79.5/mi 
(1M  cells)  and  40/un  (4M  cells).  Thus,  the  grid  spacing  along  the  surface  is  comparable 
to  the  size  of  the  heat  transfer  sensor,  but  significantly  smaller  than  the  diameter  of  the 
pressure  sensor.  The  width  of  the  symbols  used  for  the  surface  pressure  graphs  in  Section  4 
are  approximately  the  size  of  the  pressure  sensor,  and  thus  it  is  possible  to  visually  “average” 
the  computed  pressure  distribution  over  the  size  of  the  pressure  sensor. 

A  separate  extensive  set  of  experiments  were  performed  at  CUBRC  to  examine  the  shock 
wave  turbulent  boundary  layer  interactions  under  equilibrium  conditions  in  the  range  Mach  5 
to  8  for  a  cone  flare  and  hollow  cylinder  flare  of  length  2.50  m  and  2.79  m,  respectively. 
Preliminary  computations  were  performed  for  the  hollow  cylinder  flare  geometry  at  one  test 
condition;  however,  the  computational  requirements  for  the  simulation  significantly  exceeded 
the  available  resources  at  Rutgers  University,  and  therefore  efforts  were  focused  on  a  full 
study  of  the  two  laminar  shock  wave  boundary  layer  interaction  configurations  described 
above  at  four  test  conditions  each. 
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3  Methodology 


Navier-Stokes  simulations  were  performed  for  the  double  cone  and  hollow  cylinder  using  two 
different  thermochemistry  models  -  perfect  gas  and  Park  I  [12]. 

3.1  Perfect  Gas  Laminar  Navier-Stokes 

The  governing  equations  are  the  laminar  Navier-Stokes  equations  for  a  perfect  gas.  Using 
the  Einstein  summation  notation, 


dp  dpUj 
Ot  dxj 

dpiii  OpUiUj 


Ot 


dxj 


=  0 

_  _dp_  djii 

dxi  dxj 


dps  d  ,  .  dcji  drijUi 

-jr-  +  7r-(pe  +  p)uj  =  -TTL  +  ^— 
dt  dxj  dxj  dxj 

p  =  pRT 


where  the  total  energy  per  unit  mass  e  is 

e  =  e  +  \ujUj 

where  the  internal  energy  per  unit  mass  e  is 

e  =  CyT 

and  the  heat  flux  vector  and  laminar  viscous  stress  tensor  are 

,  dT 

=  -ko7t 


duh 


Tij  ~  ~^d7kSil  +  fl\~  +  ~ 


dui  dui 


dxj  dxi 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


The  molecular  viscosity  /« is  defined  by  Sutherland’s  Law  and  the  molecular  Prandtl  number 
Pr  =  pcp/k  is  0.72.  The  gas  constant  R  =  287  J/kg-K  for  air. 


3.2  Non-equilibrium  Laminar  Navier-Stokes 

We  consider  a  reacting  mixture  of  gases  with  density  pa  for  a  =  l,...,n  of  which  a  = 
1, . . . ,  m  constitute  diatomic  (or  polyatomic)  species  and  the  remainder  (i  =  m  +  1, . . . ,  n) 
represent  monatomic  species. 
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Conservation  of  Mass 
The  conservation  of  mass  is 


d£a  dpQUj  =  d_ 

dt  dxj  a  dxj 

where  pa  is  the  density  of  species  a ,  the  mass-averaged 
density 


N 


for  a  =  1, . . .  ,n  (9) 

velocity  is  Uj ,  and  p  is  the  mixture 


The  mass  fraction  is  defined  as 


P=  I] 

a=l 

(10) 

ii 

(11) 

The  rate  of  production  of  species  a  is  denoted  as  w^pe  and  defined  as 


3= 1 


iS(£rr£s(a; 


for  a  =  1, . . . ,  n  (12) 


for  the  general  reaction  expressions 


+ 


+  U'n,jXn 


+  •  •  •  + 


for  j  =  1, . . . ,  J 


(13) 


where  J  is  the  number  of  reactions,  and  u'  ,  and  v"  ■  are  the  stoichiometric  coefficients  of 
the  reactants  and  species  Xa  in  the  jth  reaction. 


The  diffusion  of  species  is  modeled  by  Fick’s  Law  assuming  a  uniform  diffusivity  D  defined 
by 


(14) 


where  Sc  =  0.7  is  the  constant  Schmidt  number  and  p  is  the  molecular  viscosity  defined  by 
Sutherland’s  Law  or  Svehla  and  McBride  [13]. 


Conservation  of  Momentum 


The  conservation  of  momentum  is 


+  o^uj  =  _dp  +  fori  =  lj23 


dt 


Oxj 


dxi  Ox, 


(15) 


where  r,j  is  the  laminar  viscous  stress  tensor  defined  by  Eq  (8). 


Conservation  of  Total  Energy 

The  total  energy  per  unit  mass  e  is  the  sum  of  the  internal  energy  per  unit  mass  e  and  the 
kinetic  energy  per  unit  mass 

e  =  e+\ujUj  (16) 
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The  internal  energy  per  unit  mass  e  is  the  sum  of  the  internal  energies  of  each  of  the  n 
species 

e=E^ea  (17) 

a=l  P 

where  the  internal  energy  per  unit  mass  of  each  species  eQ  is  the  sum  of  an  equilibrium 
internal  energy  e^q(T)  due  to  random  translational  energy  and  rotational  energy  (in  the  case 
of  molecules)  at  a  bulk  equilibrium  temperature  T  and  a  non-equilibrium  internal  energy 
e£b(TQvib)  due  to  vibrational  excitation  (in  the  case  of  molecules) 

e«  =  C(r)  +  CCCb)  (is) 

The  equilibrium  internal  energy  of  species  a  is 

e"(r)  =  +  r  cv„(T)dT  (19) 

JTtc  f 


The  conservation  of  total  energy  is 


dps  d  dTijUi  dqj 

-W  +  aT/l,e+p)u>  =  -^--dFj 


where  the  heat  transfer  vector  is  defined  by 


UXJ  a—l  ux3  a=l  ux3 


(20) 


(21) 


The  static  enthalpy  per  unit  mass  for  species  ct  is 

hn  =  h°fn  +  /  cPa  (T)  dT 

JTtc  f 

where  hja  is  the  enthalpy  of  formation  of  species  a  at  Ttc,. 
Conservation  of  Vibrational  Energy 
The  conservation  of  vibrational  energy  is 


dpe*b  dpe^uj 
dt  dxj 

The  heat  transfer  vector  is 


dqT3 

Dxj 


+  c: 


for  a  =  1 , . . . ,  m 


BTv‘b 

„vib  _  _ Z-vib  yj±  a 

Vcxj 


dxj 


-  pDae 


vib 

Q 


dYa 

dxj 


(22) 


(23) 


(24) 


The  source  term  is 


(25) 


where  e*ib  is  the  translational-vibrational  energy  transfer  per  unit  mass  of  species  a  defined 
by  the  Landau-Teller  model  [14] 


e:ib  = 


ef  (T)-<“  (T-‘) 


(20) 


where  e”ib*  is  the  equilibrium  vibrational  energy  per  unit  mass  of  species  a  defined  by 


C\T)  =  X 


R«Qtb{n) 


n=  1  exp  (0aib(7l)/T)  -  1 


(27) 


and  Are>Q  is  the  number  of  characteristic  temperatures  0^h  and  rQ  is  the  relaxation 
time  [15]  of  species  a  defined  by 


£/?= 1  Mf)Ta0 


(28) 


where  Ma  —  pa/ Ma  is  the  molar  concentration  of  species  a  and  ra0  is  the  characteristic 
relaxation  time  of  species  a  resulting  from  collisions  with  species  B  defined  by[l5] 


1  . 

Tal 3  =  ~  exp 

P 


Aa0  (r~l  -  Bn^  -  18.42 


atm-s 


where 


4 

vib  3 


AaP  =  0.00116  Ml0Qt  and  £a/?  =  0.015M^, 
and  the  averaged  molecular  weight  is  defined  by 


(29) 


(30) 


Mn0  = 


MaM0 
Ma  +  -M.fi 


Equation  of  State 
The  equation  of  state  is 

n 

V  =  T  X  PcxR* 

n= 1 

where  the  gas  constant  Rn  for  species  a  is 


Ra  — 


n 

m2 


(31) 


(32) 


(33) 


where  7 Z  is  the  Universal  Gas  Constant  and  Mn  is  the  molecular  weight  of  species  a. 
Thermochemistry  Model 

The  non-equilibrium  Navier-Stokes  simulations  used  the  thermochemistry  model  proposed 
by  Park[12]  comprising  five  species  (N2,  O2,  NO,  N,  O)  and  seventeen  reactions  (Table  3). 
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Table  3:  Thermochemistry  Model  Reactions 


Reaction 

C 

(m3/kgmolc-b) 

V 

e/k 

(k; 

a 

& 

v4i 

A2 

^3 

^4 

^5 

N2  +  N2  ->  N  +  N  +  N2 

3.70  •  101B 

-1.6 

113200 

0.5 

0.5 

10.81 

-12,61 

0  683 

-0.118 

0.006 

n2  +  n->n  +  n  +  n 

1.11  •  1019 

-1.6 

113200 

0.5 

0.5 

10.81 

-12.61 

0.683 

-0.118 

0.006 

N2  +  NO  N  +  N  +  NO 

3.70  ■  1018 

-1.6 

113200 

0.5 

0.5 

10.81 

-12.61 

0.683 

-0.118 

0.006 

N2  +  02  N  +  N  - f  02 

3.70  ■  1018 

-1.6 

113200 

0.5 

0.5 

10.81 

-12.61 

0.683 

-0.118 

0.006 

7V2  +  0->N  +  Ar  +  0 

1.11  •  1019 

-1.6 

113200 

0.5 

0.5 

10.81 

-12.61 

0.683 

-0.118 

0.006 

N2  +  0  NO  +  N 

3.18  -  IO10 

0.1 

37700 

1.0 

0  0 

2.349 

-4.828 

0.455 

-0  075 

0.004 

02  +  N2  — y  O  4-  O  4-  N2 

2.75  1016 

-1.0 

59500 

0.5 

0.5 

8.243 

-4.127 

-0.616 

0.093 

-0.005 

02  +  N  — >  O4-O4- A 

8.25  1016 

-1  0 

59500 

0.5 

0.5 

8.243 

-4.127 

-0.616 

0.093 

-0.005 

02  4-  NO  -4  0  +  0  +  NO 

2.75  1016 

-1.0 

59500 

0.5 

0.5 

8.243 

-4.127 

-0  616 

0,093 

-0.005 

02  +  02  ->  O  +  O  4-  02 

2.75  •  1016 

-1,0 

59500 

0.5 

0.5 

8.243 

-4.127 

-0.616 

0.093 

-0.005 

o2  +  0-40  +  0  +  0 

8.25  ■  1016 

-1.0 

59500 

0.5 

0.5 

8.243 

-4.127 

-0.616 

0.093 

-0.005 

NO  +  O  N  +  02 

2.16-  105 

1.29 

19220 

1.0 

0.0 

0.215 

-3.657 

0.843 

-0.136 

0.007 

NO  +  N2  ^  N  +  0  +  N2 

2.30  1014 

-0.5 

75500 

0.5 

0.5 

8.457 

-7.784 

0.228 

-0.043 

0  002 

NO  +  NO  ^  N  +  O  +  NO 

2.30-  1014 

-0.5 

75500 

0.5 

0.5 

8.457 

-7.784 

0.228 

-0.043 

0.002 

NO  4-  02  -4  N  4-  O  4-  02 

2.30  1014 

-0.5 

75500 

0.5 

0.5 

8.457 

-7.784 

0.228 

-0.043 

0  002 

N0  +  7V->7V4-0  +  7V 

4.60-  1014 

-0.5 

75500 

0.5 

0.5 

8.457 

-7.784 

0.228 

-0.043 

0.002 

NO  4-0  4  Af +  04-0 

4.60  1014 

-0.5 

75500 

0.5 

0.5 

8.457 

-7.784 

0.228 

-0.043 

0.002 

NOTES 

kf  =  CT2e-t/kT° ,  Ta  =  Ta  7*,,  ke  =  exp  (Ai  +  A2z  +  A3z2  +  A4z3  +  A 5z4)  where  2  =  10 i/T 


3.3  Numerical  Algorithm 

The  perfect  gas  laminar  Navier-Stokes  equations  were  solved  using  the  GASPex  softwafe  [16], 
and  the  11011-equilibrium  laminar  Navier-Stokes  equations  using  the  Park  I  model  were  solved 
using  the  GASP  software  [16].  The  inviscid  fluxes  were  discretized  using  either  the  algorithm 
of  Roe  [17],  Van  Leer  [18]  or  HLLE  [19]  with  Min-Mode  reconstruction.  Viscous  fluxes  were 
discretized  using  central  differencing. 

A  multi-zone  block-structured  grid  was  generated  for  both  configurations.  The  computa¬ 
tional  domains  are  shown  in  Fig.  4  and  5.  For  the  double  cone,  axisymmetry  conditions 
are  applied  from  A  to  B,  and  symmetry  conditions  on  C  to  D.  Inflow  boundary  conditions 
are  applied  from  B  to  C,  and  from  D  to  G.  No  slip  boundary  conditions  are  applied  011  H 
to  A  with  fixed  wall  temperature  and  110  surface  reactions  (for  Park  I  model  simulations1). 
Zero  gradient  boundary  conditions  are  applied  on  G  to  H.  For  the  hollow  cylinder  flare,  the 
boundary  conditions  are  extrapolation  from  A  to  B,  inflow  from  B  to  C,  symmetry  from  C 
to  D.  inflow  form  D  to  E,  symmetry  from  E  to  F,  inflow  from  F  to  G,  symmetry  from  G  to 
H,  extrapolation  from  II  to  I,  isothermal  no-slip  wall  from  I  to  J  and  then  to  A. 

For  the  double  cone,  typically  two  computations  were  performed  for  each  configuration  and 
inflow  condition.  The  first  grid  was  typically  1  M  cells,  and  the  second  grid  was  4  M  or  5.3  M 
cells.  For  the  hollow  cylinder  flare,  three  computations  were  performed  for  each  perfect  gas 
simulation  (except  for  the  21.85  M.T/kg  case)  at  1  M,  4.2  M  and  16.7  M  cells.  Two  simulations 
were  performed  using  the  Park  I  for  the  10.43  MJ/kg,  15.54  MJ/kg  and  21.85  MJ/kg  cases 

1 A  separate  simulation  was  performed  using  air  catalysis  boundary  conditions  on  the  boundary  H  to  A 
and  the  results  were  essentially  identical  to  the  no  surface  reaction  case. 
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with  1  M  and  4  M  cells. 

The  simulations  were  initialized  using  the  freestream  conditions  and  converged  to  steady  state 
using  Gauss-Seidel  iteration  or  Dual  Time  Stepping  [16].  Computations  were  performed  on 
48  or  128  core  Linux  clusters  in  the  Department  of  Mechanical  and  Aerospace  Engineering, 
and  on  112  core  Linux  cluster  in  the  School  of  Engineering. 


F  G 


H 


axis 


Figure  4:  Double  Cone 


G  H 


Figure  5:  Hollow  Cylinder  Flare 
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4  Results 


An  extensive  series  of  simulations  were  performed  for  the  double  cone  and  hollow  cylinder 
flare  geometries  using  both  the  perfect  gas  and  Park  I  thermochemisty  models.  Addition¬ 
ally,  a  grid  refinement  study  was  performed  for  selected  cases  to  assess  the  accuracy  of  the 
simulations.  The  results  are  presented  below. 

4.1  Double  Cone 


Perfect  gas  and  Park  I  thermochemistry  model  simulations  were  performed  for  all  four  ex¬ 
perimental  test  conditions  listed  in  Table  1.  Details  are  presented  in  Table  4. 


Table  4:  Double  Cone  Simulations 


Run  Total  Mach 

No  Enthalpy  Number 

( MJ/kg) 

1  5.44  12.2 

2  9.65  10.9 

6  15.23  11.46 

4  21.77  12.82 


Perfect  Park  I 
Gas 

1M  4M  1.3M  5.3M 


V  V  V 

V  V 

V  V  V 

V  si 


The  basic  flowfield  structure  for  21.77  MJ/kg  (Run  4)  is  displayed  in  Fig.  6  for  reference 
in  comparison  of  the  experiment  and  simulations.  A  detailed  discussion  of  the  flowfield 
structure  is  presented  later. 


(a)  Double  cone 


(b)  Corner  region 


Figure  6:  Mach  contours  for  double  cone  using  Park  I  (21.77  MJ/kg) 
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4.1.1  5.44  MJ/kg  (Run  1) 


The  computed  and  experimental  heat  transfer  and  surface  pressure  for  5.44  MJ/kg  (Run  1) 
using  the  perfect  gas  and  Park  I  models  are  presented  in  Fig.  8.  The  perfect  gas  and 
Park  I  results  are  in  generally  close  agreement  with  each  other  for  both  surface  pressure 
and  heat  transfer  except  for  the  location  of  the  peak  heat  transfer  and  peak  pressure.  The 
computed  heat  transfer  upstream  of  the  separation  point  at  x  =  7  cm  slightly  underpredicts 
the  experiment.  The  location  of  separation  as  defined  by  the  drop  in  Qw  is  accurately 
predicted  by  both  models. 

The  computed  peak  heat  transfer  and  pressure  overestimate  the  experiment  by  40%  and 
60%,  respectively,  for  the  perfect  gas  model.  The  computed  peak  heat  transfer  for  the  Park  I 
model  is  within  the  experimental  uncertainty;  however,  the  computed  peak  pressure  for  the 
Park  I  model  is  55%  above  the  experiment. 

A  grid  refinement  study  was  performed  for  5.44  MJ/kg  (Run  1)  using  the  perfect  gas  model. 
Results  for  the  1  M  and  4  M  cell  grids  (Fig.  8)  are  virtually  identical,  with  the  maximum 
difference  in  peak  heat  transfer  and  peak  pressure  between  the  two  grid  solutions  of  6%  and 
3%,  respectively. 

4.1.2  9.65  MJ/kg  (Run  2) 

The  computed  and  experimental  heat  transfer 
and  surface  pressure  for  9.65  MJ/kg  (Run  2) 
using  the  perfect  gas  and  Park  I  models  are 
presented  in  Fig.  9.  The  perfect  gas  results 
show  significant  disagreement  with  the  experi¬ 
ment.  The  streamwise  length  of  the  separated 
region  is  overpredicted  by  75%.  The  mag¬ 
nitude  of  the  peak  heat  transfer  is  in  good 
agreement  with  experimental  peak  heat  trans¬ 
fer;  however,  the  location  is  shifted  down¬ 
stream  by  1.5  cm.  The  Park  I  results  show 
closer  agreement  with  experiment  than  the 
perfect  gas  simulations,  although  significant 
differences  are  evident.  The  streamwise  length 
of  the  separation  region  is  underpredicted  by 
40%.  The  peak  heat  transfer  is  predicted 
within  15%,  while  the  computed  peak  pres¬ 
sure  is  30%  greater  than  the  peak  experimen¬ 
tal  pressure. 


Figure  7:  Temperature  contours  (Park  I) 
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The  dramatic  drop  and  rise  in  the  computed  heat  transfer  at  x  =  10.05  cm  using  Park  I  is  a 
real  physical  phenomenon  in  the  simulation.  The  numerical  grid  spacing  along  the  surface  in 
this  region  (see  inset  in  Fig.  9(d))  indicates  a  careful  resolution  of  the  heat  transfer  profile. 
This  phenomenon  is  associated  with  the  formation  of  a  small  recirculation  bubble  on  the 
wall  as  seen  in  Fig.  7. 

4.1.3  15.23  MJ/kg  (Run  6) 

The  computed  and  experimental  heat  transfer  and  surface  pressure  for  15.23  MJ/kg  (Run  G) 
using  the  perfect  gas  and  Park  I  models  are  presented  in  Fig.  10.  Similar  to  the  9.G5  MJ/kg 
case,  the  perfect  gas  results  show  significant  disagreement  with  the  experiment.  The  stream- 
wise  length  of  the  separation  region  and  the  peak  heat  transfer  are  both  overpredicted  by 
100%. 

The  Park  I  results  show  closer  agreement  with  experiment  than  the  perfect  gas  simulations, 
although  significant  differences  are  evident  as  seen  previously  for  the  9.G5  MJ/kg  case.  The 
streamwise  length  of  the  separation  region  is  underpredicted  by  50%.  The  magnitude  of  the 
peak  heat  transfer  is  accurately  predicted. 

A  grid  refinement  study  was  performed  for  15.23  MJ/kg  (Run  6)  using  the  Park  I  model. 
Results  for  the  1.3M  and  5.3M  cell  cases  show  very  close  agreement  (Fig.  10). 

The  Park  I  heat  transfer  and  pressure  profiles  display  a  dramatic  drop  and  rise  at  x  = 
10.05  cm  similar  to  the  9.65  MJ /kg  case.  The  phenomenon  is  captured  by  both  the  1.3M  and 
5.3M  cell  grid  systems.  The  behavior  is  attributable  to  the  formation  of  a  small  recirculation 
bubble  on  the  surface  near  reattachment. 

4.1.4  21.77  MJ/kg  (Run  4) 

The  computed  and  experimental  heat  transfer  and  surface  pressure  for  21.77  MJ/kg  using 
the  perfect  gas  and  Park  I  models  are  presented  in  Fig.  11.  Similar  to  the  9.65  MJ/kg  and 
15.23  MJ /kg  cases,  the  perfect  gas  results  show  significant  disagreement  with  the  experiment. 
The  streamwise  length  of  the  separation  region  and  the  magnitude  of  the  peak  heat  transfer 
are  overpredicted  by  100%.  The  Park  I  results  show  closer  agreement  with  experiment  than 
the  perfect  gas  simulations,  although  minor  differences  are  evident.  The  streamwise  length 
of  the  separation  region  is  accurately  predicted,  together  with  the  magnitude  of  the  peak 
heat  transfer  and  peak  pressure. 
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(W/cm‘)  p  (kPa) 


(a)  Pressure  (perfect  gas)  (b)  Pressure  (Park  I) 


(c)  Heat  Transfer  (perfect  gas) 


(d)  Heat  Transfer  (Park  I) 


Figure  8:  Surface  Pressure  and  Heat  Transfer  for  Run  1  (5.44  M J /kg) 
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Q*  (W/cm*)  p  (kPa) 


(a)  Pressure  (perfect  gas) 


(c)  Heat  Transfer  (perfect  gas) 


(b)  Pressure  (Park  I) 


(d)  Heat  Transfer  (Park  I) 


Figure  9:  Surface  Pressure  and  Heat  Transfer  for  Run  2  (9.65  MJ/kg) 
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Q*  (VWcirn  ^  w  p  (kPa) 


(a)  Pressure  (perfect  gas)  (b)  Pressure  (Park  I) 


(c)  Heat  Transfer  (perfect  gas) 


(d)  Heat  Transfer  (Park  I) 


Figure  10:  Surface  Pressure  and  Heat  Transfer  for  Run  6  (15.23  MJ/kg) 
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(a)  Pressure  (perfect  gas)  (b)  Pressure  (Park  I) 


(c)  Heat  Transfer  (perfect  gas)  (d)  Heat  Transfer  (Park  I) 


Figure  11:  Surface  Pressure  and  Heat  Transfer  for  Run  4  (21.77  MJ/kg) 
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4.1.5  Analysis  of  Flowfield 


Fig.  12  displays  contours  of  species  mass  fractions  of  NO,  O  and  N  for  the  9.65  MJ/kg  case2. 
A  significant  fraction  of  atomic  oxygen  O  is  created  by  the  aft  cone  shock  as  indicated  in 
Figs.  12(e,f).  The  maximum  value  Y0  —  0.157,  and  the  maximum  value  Y0/Y0i  =  4.67.  The 
mass  fraction  of  atomic  nitrogen  N  is  very  low,  with  a  maximum  mass  fraction  yN  =  0.006, 
and  hence  is  not  visible  in  Figs.  12(c.d).  The  atomic  nitrogen  N  generated  by  dissociation 
of  N2  at  the  aft  shock  rapidly  reacts.  The  maximum  mass  fraction  TNO  =  0.084. 

Fig.  13  displays  contours  of  species  mass  fractions  of  NO,  O  and  N  for  the  15.23  MJ/kg 
case.  Similar  to  the  9.65  MJ/kg  case,  a  significant  fraction  of  atomic  oxygen  O  is  created  by 
the  aft  cone  shock  as  indicated  in  Figs.  13(e,f).  The  maximum  value  Y0  =  0.2346,  and  the 
maximum  value  Y0/Y02  =  1942.  Nearly  all  of  the  molecular  oxygen  O2  is  dissociated  by  the 
high  temperatures  generated  by  the  aft  cone  shock,  with  the  mass  fraction  YC)2  substantially 
below  1%  downstream  of  the  aft  cone  shock.  The  mass  fraction  of  atomic  nitrogen  N  increases 
downstream  of  the  aft  cone  shock  in  comparison  with  the  9.65  M J /kg  case,  with  a  maximum 
mass  fraction  TN  =  0.084.  The  maximum  mass  fraction  TNO  =  0.088. 

Fig.  14  displays  contours  of  species  mass  fractions  of  NO,  O  and  N  for  the  21.77  MJ/kg 
case.  Similar  to  the  9.65  MJ/kg  and  15.23  MJ/kg  cases,  a  significant  fraction  of  atomic 
oxygen  O  is  created  by  the  aft  cone  shock  as  indicated  in  Figs.  14(e,f).  The  maximum  value 
Y0  —  0.2349,  and  the  maximum  value  Y0/Y02  =  2.60- 104,  indicating  effectively  all  molecular 
oxygen  O2  has  been  dissociated  at  some  point  downstream  of  the  aft  cone  shock.  The  mass 
fraction  of  atomic  nitrogen  N  indcreases  downstream  of  the  aft  cone  shock  in  comparison 
with  the  9.65  MJ/kg  and  15.23  MJ/kg  cases,  with  a  maximum  mass  fraction  TN  =  0.184. 
The  maximum  mass  fraction  TNO  =  0.093. 

The  difference  between  the  vibrational  temperature3 and  translational-rotational  temperature 
Tvib— T  for  the  9.65  MJ/kg,  15.23  MJ/kg  and  21.77  MJ/kg  cases  is  displayed  in  Fig.  15.  In  the 
incoming  flow,  Tvib  =  T.  The  vibrational  temperature  remains  relatively  frozen  throughout 
the  interaction,  and  consequently  Tvib  —  T  is  negative  immediately  downstream  of  the  aft 
cone  shock,  and  positive  downstream  of  the  expansion  corner. 


2The  mass  fraction  scale  is  the  same  for  all  cases  and  species. 

3The  averaged  vibrational  temperature  7’"^,’  is  determined  assuming  an  equilibrium  vibrational  distribu¬ 
tion  at  the  averaged  vibrational  energy  defined  as 

m 

C  =  Ey«e«b 

Q  =  1 
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(b)  Fno  and  Mach  isolines 


(d)  Fn  and  Mach  isolines 


Figure  12:  Species  mass  fractions  for  Park  I  (9.65  MJ/kg) 


23 


(a)  Yno 


(b)  yNO  and  Mach  isolines 


(c)yN 


(d)  yN  and  Mach  isolines 


Figure  13:  Species  mass  fractions  for  Park  I  (15.23  MJ/kg) 
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(c)  rN 


(d)  Yn  and  Mach  isolines 


(f)  Y0  and  Mach  isolines 


Figure  14:  Species  mass  fractions  for  Park  I  (21.77  M.T/kg) 
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(a)  9.65  MJ/kg 


(b)  9.65  MJ/kg 


(c)  15.23  MJ/kg 


(d)  15.23  MJ/kg 


(e)  21.77  MJ/kg 


(f)  21.77  MJ/kg 


Figure  15:  Tvib  —  T  (left)  and  Tvib  —  T  with  Mach  isolines  (right) 
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4.2  Hollow  Cylinder  Flare 


Table  5:  Hollow  Cylinder  Flare  Simulations 


Run 

No 

Total 

Enthalpy 

(MJ/kg) 

Mach 

Number 

1M 

Perfect 

Gas 

4M 

16M 

Park  I 

1M  4M 

1 

5.07 

11.3 

% 

% 

V 

2 

10.43 

12.6 

% 

% 

% 

% 

4 

15.54 

11.5 

% 

% 

V 

V 

5 

21.85 

13.2 

% 

% 

% 

4.2.1  5.07  MJ/kg  (Run  1) 

The  computed  and  experimental  surface  pressure  and  heat  transfer  for  Run  1  (5.07  MJ/kg) 
are  presented  in  Fig.  16.  Figs.  16(a)  and  16(b)  display  the  results  using  the  perfect  gas  model 
including  grid  refinement.  The  computed  peak  pressure  using  the  perfect  gas  model  exceeds 
the  experiment  by  60%.  The  computed  peak  heat  transfer  using  the  perfect  gas  model  agrees 
with  experiment  within  the  experimental  uncertainty;  however,  the  perfect  gas  simulations 
overpredict  the  size  of  the  separation  region.  The  perfect  gas  simulations  demonstrated  grid 
independence  as  indicated  by  the  results  for  the  three  different  grids. 

4.2.2  10.43  MJ/kg  (Run  2) 

The  computed  and  experimental  surface  pressure  and  heat  transfer  for  Run  2  (10.43  MJ/kg) 
are  presented  in  Fig.  17.  Figs.  17(a)  and  17(c)  display  the  results  using  the  perfect  gas  model 
including  grid  refinement.  Figs.  17(b)  and  17(d)  display  the  results  of  non-equilibrium  gas 
using  the  Park  I  model  for  1M  cells.  The  computed  peak  pressure  using  the  perfect  gas  model 
exceeds  the  experiment  by  30%,  while  the  computed  peak  pressure  of  the  Park  I  exceeds  the 
experiment  by  50%.  The  computed  peak  heat  transfer  using  the  perfect  gas  model  agrees 
with  experiment  within  the  experimental  uncertainty,  while  the  calculated  peak  heat  transfer 
of  Park  I  model  exceeds  the  experiment  by  25%.  No  separated  region  is  observed  in  any  of 
the  simulations  of  the  perfect  gas  model  which  is  the  same  as  what  is  seen  in  the  experiment. 
However,  the  simulations  with  Park  I  model  creates  a  separation  region. 

4.2.3  15.54  MJ/kg  (Run  4) 

The  computed  and  experimental  surface  pressure  and  heat  transfer  for  Run  4  (15.54  MJ/kg) 
are  presented  in  Fig.  18.  Figs.  18(a)  and  18(c)  display  the  results  using  the  perfect  gas  model 
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including  grid  refinement  for  three  sequence  of  grids.  Figs.  18(b)  and  18(d)  display  the  results 
of  non-equilibrium  gas  using  the  Park  I  model  including  grid  refinement  for  two  sequence 
of  grids.  The  computed  peak  pressure  using  both  the  perfect  gas  model  and  Park  I  model 
exceeds  the  experiment  by  23%.  The  computed  peak  heat  transfer  using  the  perfect  gas 
model  agrees  with  experiment  within  the  experimental  uncertainty,  while  the  Park  I  model 
simulation  overestimate  the  peak  heat  transfer  by  27%.  Moreover,  the  separated  region  is 
overpredicted  by  the  perfect  gas  model  while  it  underpredicted  by  the  Park  I  model. 

4.2.4  21.85  MJ/kg  (Run  5) 

The  computed  and  experimental  surface  pressure  and  heat  transfer  for  Run  5  (21.85  MJ/kg) 
are  presented  in  Fig.  19.  Figs.  19(a)  and  19(c)  display  the  results  using  the  perfect  gas  model 
including  grid  refinement  for  two  sequences  of  grids,  and  Figs.  19(b)  and  19(d)  display  the 
results  using  the  non-equilibrium  model  with  Park  I  reactions.  The  computed  peak  pres¬ 
sure  using  the  perfect  gas  model  exceeds  the  experiment  by  87%,  while  the  computed  peak 
pressure  using  the  non-equilibrium  model  exceeds  the  experiment  by  47%.  The  computed 
peak  heat  transfer  using  the  perfect  gas  model  is  within  25%  of  the  experiment  while  the 
computed  peak  heat  transfer  using  the  non-equilibrium  model  agrees  with  experiment  within 
the  experimental  uncertainty.  The  separated  region  in  the  simulations  using  the  perfect  gas 
model  is  larger  than  experimental  data;  however  the  separated  region  using  the  Park  I  model 
is  in  good  agreement  with  experiment. 


(a)  Pressure  (perfect  gas) 


(b)  Heat  transfer  (perfect  gas) 


Figure  1C:  Surface  pressure  and  heat  transfer  for  Run  1  (5.07  MJ/kg) 
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(ed»)  d 


(a)  Pressure  (perfect  gas) 


(b)  Pressure  (Park  I) 


(c)  Heat  transfer  (perfect  gas) 


Figure  17:  Surface  pressure  and  heat  transfer  for  Run  2  (10.43  MJ/kg) 
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(W/cm2)  p  (kPa) 


O  Experiment 


(a)  Pressure  (perfect  gas) 


O  Experiment 


(b)  Pressure  (Park  I) 


O  Experiment 

-  Calculated  (Perfect  Gas,  1M  cells) 

-  Calculated  (Perfect  Gas,  4M  cells) 

-  Calculated  (Perfect  Gas,  16M  cells) 


(c)  Heat  transfer  (perfect  gas) 


(d)  Heat  transfer  (ParkI) 


Figure  18:  Surface  pressure  and  heat  transfer  for  Run  4  (15.54  MJ/kg) 
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Qw  (W/cm^)  p  (kPa> 


(a)  Pressure  (perfect  gas)  (b)  Pressure  (Park  I) 


(c)  Heat  tranfer  (perfect  gas) 


(d)  Heat  Transfer  (Park  I) 


Figure  19:  Surface  pressure  and  heat  transfer  for  Run  5  (21.85  MJ/kg) 
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4.2.5  Analysis  of  Flowfield 


Fig.  20  displays  the  Mach  number  contours  and  streamlines  for  Run  1  (5.07  MJ/kg).  Since 
the  5.07  MJ/kg  enthalpy  is  likely  less  than  what  is  needed  for  the  gas  to  dissociate,  only  the 
perfect  gas  model  is  used  for  simulation.  A  recirculation  region  forms  at  the  cylinder-flare 
juncture  due  to  the  adverse  pressure  gradient.  The  boundary  layer  displacement  shock,  flare 
shock  and  separation  shock  intersect  at  approximately  the  location  of  reattachment . 

Fig.  21  shows  the  Mach  number  contours  and  streamlines  for  Run  2  (10.43  MJ/kg).  A  sur¬ 
prising  result  is  the  absence  of  a  recirculation  region  at  the  corner  for  perfect  gas  model  while 
the  recirculation  region  exist  in  non-equilibrium  model.  In  the  non-equilibrium  simulation 
the  boundary  layer  displacement  shock,  flare  shock  and  separation  shock  intersect  at  about 
reattachment  point. 

Fig.  22  displays  the  Mach  number  contours  and  streamlines  for  Run  4  (15.54  MJ/kg).  A 
separated  region  is  evident  and  is  significantly  larger  than  for  Run  1  (5.07  MJ/kg)  in  the  case 
of  using  the  perfect  gas  model,  moreover,  recirculation  region  is  much  more  smaller  in  the 
non-equilibrium  model.  For  the  perfect  gas  model,  the  separation  shock  and  boundary  layer 
displacement  shock  intersect  over  the  recirculation  region  and  the  resulting  shock  intersects 
the  Flare  shock  in  the  region  above  the  reattachment  point.  However,  in  the  non-equilibrium 
simulation  the  boundary  layer  displacement  shock  and  separation  shock  don’t  intersect  with 
each  other  at  all  and  they  intersect  with  flare  shock  in  different  location. 

Fig.  23  shows  the  Mach  number  contours  and  streamlines  for  Run  5  (21.85  MJ/kg).  The 
separation  region  is  reduced  in  size  compared  to  Run  4  in  the  perfect  gas  model  while 
the  separation  region  is  completely  vanished  in  non-equilibrium  simulation  with  the  Park  I 


model. 


Mach-Number:  1  2345  678  9  10  11 


(a)  Mach  contours  (perfect  gas) 


(b)  Corner  region  (perfect  gas) 


Figure  20:  Flow  Structure  of  Run  1  (5.07  MJ/kg) 
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Mach-Number:  1  23456789101112 


Mach-Number;  1  2  3  4  5  6  7  8  9101112 


(a)  Mach  contours  (perfect  gas) 


(b)  Mach  contours  (Park  I) 


Mach-Number:  1  2 


Separation  Point 
Boundary  Layer 


Boundary  Layer 


Point 


(c)  Corner  region  (perfect  gas) 


(d)  Corner  region  (Park  I) 


Figure  21:  Flow  Structure  of  Run  2  (10.43  MJ/kg) 


33 


(a)  Mach  contours  (perfect  gas) 


(b)  Mach  contours  (Park  1) 


Reattachment 
Point 

Recirculation 

Separation  Point 
Boundary  Layer 


(c)  Corner  region  (perfect  gas) 


(d)  Corner  region  (Park  I) 


Figure  22:  Flow  Structure  of  Run  4  (15.54  MJ/kg) 
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Mach-Number:  1  23456789  10  11  12 


Mach-Number:  1  23456789101112 


(a)  Mach  contours  (perfect  gas) 


(b)  Mach  contours  (Park  1) 


Mach-Number  1  2  3  4  5  6  7  8  9  10  11  12 


Mach-Number:  1  2  3  4  5  6  7  6  9  10  11  12 


Reattachment 
Point 

Recirculation 


Separation  Point 
Boundary  Layer 


Boundary  Layer 


(c)  Corner  region  (perfect  gas) 


(d)  Corner  region  (Park  I) 


Figure  23:  Flow  Structure  of  Run  5  (21.85  MJ/kg) 

Contours  of  mass  fractions  for  N,  O  and  NO  are  displayed  in  Fig.  24  to  26.  Fig.  24  displays 
the  dissociation  of  N2  and  02  occur  in  the  junction  of  cylinder  and  flare.  The  mass  fraction 
of  N  is  very  small  in  relative  to  mass  fraction  of  NO  and  O  and  occurs  mostly  along  the  flare. 
Fig.  25  demonstrates  the  significant  increase  in  the  dissociation  of  N2  and  02  molecules  by 
increasing  total  enthalpy  from  10.43  MJ/kg  in  Run  2  to  15.54  MJ/kg  in  Run  4.  In  this 
case,  most  of  the  mass  fraction  of  NO  molecule  and  N  and  O  atoms  appear  along  the  flare 
especially  near  the  reattachment  point.  Fig.  26  displays  significant  dissociation  of  N2  and 
02  occur  downstream  of  the  reattachment  point,  and  NO  rapidly  forms.  Figs.  28  and  29 
display  Tvih  —  T.  The  vibrational  temperature  significantly  lags  the  translational-rotational 
temperature  downstream  of  the  reattachment  point. 
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(a)  Yno 


(b)  yNO  with  constant  Macli  lines 


Figure  24:  Mass  fraction  of  NO,  N  and  O  contours  for  Run  2  (10.43  MJ/kg) 
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(e)  Y0 


(f)  Y0  with  constant  Mach  lines 


Figure  25:  Mass  fraction  of  NO,  N  and  O  contours  for  Run  4  (15.54  MJ/kg) 
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(c)  yN 


(e)YQ 


(d)  yN  with  constant  Mach  lines 


(f)  Y0  with  constant  Mach  lines 


Figure  26:  Mass  fraction  of  NO,  N  and  0  contours  for  Run  5  (21.85  MJ/kg) 
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(a)  Tvib  —  T  contour  (b)  Tvlb  —  T  with  constant  Mach 

lines 

Figure  27:  Vibrational  temperature  contours  with  Mach  lines  for  Run  2  (10.43  MJ/kg) 


Tv-T:  -5470  -4009  -2549  -1088  372 


Tv-T:  -5470  -MW  -2549  1068  372 


lines 


Figure  28:  Vibrational  temperature  contours  with  Mach  lines  for  Run  4  (15.54  MJ/kg) 


Tv  T:  -7079  -5229  -3378  -1528  322 


Tv-T:  -7079  -6229  -3378  -1528  3 72 


Figure  29:  Vibrational  temperature  contours  with  Mach  lines  for  Run  5  (21.85  MJ/kg) 
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5  Conclusions,  Significance  of  Results  and  Transitions 

The  principal  conclusions  of  the  research  are 

•  The  Park  I  model  accurately  predicts  the  peak  heat  flux  max  Qw  in  nearly  all  cases 

•  The  Park  I  model  accurately  predicts  the  separation  length  Lecp  and  peak  pressure 
max  p  in  approximately  one-half  of  the  cases 

•  The  perfect  gas  model  accurately  predicts  the  peak  heat  flux  max  Qw  in  all  four  hollow 
cylinder  flare  cases,  but  in  none  of  the  double  cone  cases 

•  The  perfect  gas  model  accurately  predicts  the  separation  length  Laop  in  one  case  and 
peak  pressure  max  p  in  none  of  the  cases 


Table  6:  Accuracy  of  Perfect  Gas  and  Park  I  Models 


Geometry 

Run 

No. 

Total 

Enthalpy 

(MJ/kg) 

Perfect  Gas 
max  p  max  Qw 

-^sop 

Park  I 
max  p 

max  Qw 

1 

5.44 

• 

o 

o 

• 

o 

• 

Double 

2 

9.65 

o 

o 

o 

o 

o 

• 

Cone 

6 

15.23 

o 

o 

o 

o 

• 

• 

4 

21.77 

o 

o 

o 

• 

• 

• 

Hollow 

1 

5.C7 

o 

o 

• 

n/a 

n/a 

n/a 

Cylinder 

2 

10.43 

o 

o 

• 

• 

O 

• 

Flare 

4 

15.54 

o 

o 

• 

• 

• 

• 

6 

21.85 

o 

o 

• 

• 

O 

• 

NOTES: 

•  Accurate 
O  Inaccurate 

Lscp  is  the  streanrwise  separation  length  defined  by  the  Qw  profile 
max  p  and  max  Qw  refer  to  the  peak  values  of  pressure  and  heat  transfer 


The  significance  of  the  results  is: 

•  The  Park  I  model  is  a  reliable  model  for  prediction  of  peak  heat  flux 

•  The  Park  I  model  is  not  a  reliable  model  for  prediction  of  peak  surface  pressure  or 
separation  length  ( i.e .,  surface  pressure  distribution) 


40 


•  The  perfect  gas  model  is  not  a  reliable  model  for  prediction  of  peak  surface  pressure  or 
separation  length  (z.e.,  surface  pressure  distribution)  except  at  the  lowest  stagnation 
enthalpy 

•  The  perfect  gas  model  is  not  a  reliable  model  for  prediction  of  peak  heat  flux  except  at 
the  lowest  stagnation  enthalpy;  it  is  accurate  for  hollow  cylinder  flare  at  all  stagnation 
enthalpies,  but  not  for  the  double  cone 

The  significance  is  limited  to  the  types  of  shock  wave  laminar  boundary  layer  interactions 
considered  in  this  study. 
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